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Decoherence free subspaces (DFS) is a theoretical tool towards experimental implementation of 
quantum information storage and processing. However, they represent an experimental challenge, 
since conditions for their existence are very stringent. This work explores the situation in which a 
system of N oscillators coupled to a bath of harmonic oscillators is close to satisfy the conditions for 
the existence of DFS. We show, in the Born-Markov limit and for small deviations from separability 
and degeneracy conditions, that there are weak decoherence subspaces which resemble the original 
notion of DFS. 

PACS numbers: 03.65.Yz,03.67.Lx 

I. INTRODUCTION 

The very same mechanism responsible for the potential improvements on computation speed using quantum mechan- 
ics, is the one which greatly hinders immediate technical implementation. Entanglement between different subsystems 
is essential for the production of the states used in information processing [|; at the same time, as these qubits can 
not be completely isolated from its environment, entanglement with the environmental degrees of freedom is a gen- 
eral feature. The deleterious effect of this coupling is usually called decoherence)^. Therefore much effort has been 
devoted to finding ways around decoherence in quantum computation, such as error correcting codesQ, dynamical 
decoupling^ and computation in decoherence free subspaces 0, 0- Experimental observations of decoherence free 
evolution have been reportedjl, ||. Many physical implementations have been proposed including cavity QEP[Io|. 
ions traps [lj, nuclear magnetic resonance |l2fl and semiconductor quantum dots|l3]. From the theoretical point of 
view, recent work has been mainly on proving the existence of DF subspaces, in general related to symmetries of 
the system which are preserved by the interaction with the environment, on searching for mechanisms of dynamical 
creation of DFS [3, and on the analysis of their robustness[l5j. More realistic modelsyj are scarce, and fail to provide 
insight on the effects of slightly relaxing the conditions necessary for the existence of DFS. 

In the present work we consider the case of N independent oscillators linearly coupled to a single environment, 
and show that the existence of strict decoherence free subspaces can be obtained under the following two conditions: 
degeneracy of the oscillators, and separability of the coupling with the environment. Both can be viewed as conse- 
quences of symmetries: the first involving only the system, and the second the interaction. The exact form of the 
spectral density and the temperature of the environment are immaterial in what concerns the existence of DFS, since 
they really decouples from the environment. Master equations for the evolution of the reduced density matrix of 
the system are derived with and without the Born-Markov approximation, and only the coefficients vary from one 
case to the other. For two independent oscillators we solve the dynamics of the reduced density and explicit the 
decoherence free subspace. Also in the case of two harmonic oscillators we study the effect of relaxing the degeneracy 
and separability conditions. We verify that there is no more DFS, but there remains a long leaved mode, which we 
call weak decoherence mode, and its counterbalance, a strong decoherence mode also appear. The time scales for the 
duration of these components are derived in terms of the appropriate parameters. We analyze these findings in the 
context of the robustness proof for DFS presented in Ref . . 

This contribution is organized as follows: in section[n]an introduction to the concept of decoherence free subspaces 
is given, in section ITTT1 the model with oscillators is described, and the decoherence free modes are exhibited. A short 
discussion of the classical counterpart of decoherence free modes is made. Section IIVI is devoted to the derivation of 
master equations. In section we discuss in details the case of two harmonic oscillators. As a simplifying tool, we 
introduce the notion of superoperator. The section IVT1 discuss the case of small departures from the degeneracy and 
separability conditions, and shows that despite the fact that the concept of DFS is no more applicable, there remains 
a weak decoherence mode which can be useful for quantum information storing. We give some concluding remarks. 
Some intermediate calculations have been relegated to the appendix. 
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II. A QUICK WAY TO DECOHERENCE FREE SUBSPACES 

The notion of decoherence free subspaces (DFS) can be easily captured by considering a special coupling between 
a system and its environment. Consider a system with its autonomous Hamiltonian Hg, an environment described 
by Hb, and the interaction between them given by H/. The complete Hamiltonian is then 

H = H S + H £ + H 7 . (1) 

Suppose the interaction term can be written in the form H/ = A.s'Be, with As (resp. Bg) acting only on the system 
(resp. environment) degrees of freedom. We will call this a separability condition. Any common eigenvector \a) of As 
and Hg (with eigenvalues a and h a ) does not get entangled with the environment, since in this case 



H 



| a) ® |e) = (H s |a)) ® |e) + |a) ® (ft £ |e>) + (A s |a)) ® (b e |e)) 
= |a) (h a + ft E +aB E ) \e) . 



(2) 



By linearity, any common eigenspace of As and Hs is a DFS of the system. 

Degeneracy is generally originated by symmetry. Therefore, if one finds a symmetric system where interaction with 
the environment preserves this symmetry, then any eigenspace of the system is a DFS. 

In the language of Ref . [£j , As is the only error generator, and (common) eigenspaces of (all) error generators are 
DFS. An important distinction that we make is to consider as DFS only the common eigenspaces of As and Hs, i.e.: 
the system Hamiltonian should not take the state out of a DFS. 

III. A MODEL WITH OSCILLATORS 

We now present a different situation in which DFS can be achieved. The system will consist of N identical harmonic 
oscillators (frequency u>, annihilation operators a,, we use h = 1). The environment will be modelled as a huge set of 
harmonic oscillators (frequencies cok, annihilation operators bfc). Linear coupling will be considered and the rotating 
wave approximation applied. This model is both: simple enough to be studied in details and general enough to keep 
the characteristic behaviour of the problem. It is also adequated to make the link with experimental implementations: 
one can consider vibronic states of N ions trapped together, a system in which (approximate) DFS has already been 
demonstrated |oL ITt^ . or modes of distinct cavities [l7j. or even, for N = 2, two degenerated modes of one cavity. The 
Hamiltonian to be considered is 

N 

H = ojY^ + E ^khlbk + ^(c^ajbfc + fltaa<b£). (3) 

i— 1 k i,k 

As in the previous model, we need an additional assumption on the form of the interaction. Assume the coupling 
constants gtk can be factorized as GiDk- This can be interpreted as supposing that all oscillators feel the environment 
in the same way, possibly with only a difference in strength, which depends only on the oscillator itself, not on the 
environment (the most usual models consider Gi — G, which is a special form of the here proposed model). With this 
factorizability hypothesis, one can rewrite the Hamiltonian JSJ as 



h = w J2 + E + E ( D ^E G * a i)^ + ^(E J 

i—l k k \ i i / 



(4) 



By defining the collective operators 



it takes the form 

H = ^E + E Wfc kib fe + A\ J2 c k h k + Ai 4bi, (6) 



N 
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where it is clear that only the collective mode Ai is coupled to the environment (in the above formula Ck = 
J2i \Gi\ 2 D^). One can consider the mode A x as the first of a new set of normal modes {-A-j j, an d the remain- 
ing modes thus constitute an infinite dimensional DFS. With this new set of variables, the Hamiltonian is finally 
written as 

N 

h = ^ E M&i + E + At i E c ^ + A i E c *M- ( ? ) 

i—l h k k 

While the situation in the previous section is completely quantum mechanical, the present model does have a 
classical analog, because the manipulation above can also be done with classical oscillators. In fact, there is a very 
old classical situation to which this analysis can be applied: synchronization of pendular clocks. It is known that two 
clocks in the same wall tend to synchronize in anti-phase. Each clock can be considered as an oscillator, and their 
coupling to the environment can be considered in terms of the two normal modes (phase and anti-phase modes). The 
in phase mode couples (much more) strongly to the environment, and this causes the synchronization. 

Another consistent analogy of the above model is with the superradiance in the Dicke model|l8|. In this case, 
N two-level atoms are coupled to one field mode. In the regime in which the atoms are collectively coupled to the 
field (i.e.: the field does not distinguish which atom has emitted), the radiative process can be stronger compared to 
individual emissions (due to interference). The counterpart of this process is the subradiance, i.e.: other collective 
states with weaker emission than the individual contributions (destructive interference). DFS can thus be compared 
to subradiant states of the system. 



IV. MASTER EQUATION 



As one usually does not have control on the environment degrees of freedom, the natural approach to this problem 
is to study the reduced dynamics of the N oscillators. A long but straightforward procedure 19] can be applied to 
derive the master equation 



dp _ 1 
dt ih 



H ,/5 



(A + e) ( 2A 1/ 3A T 1 - AjAxp - pA{A ± 



e ( 2AJ/3A1 - AiAjp - pAiAj 



(8) 



where 



N 



H = fi^Y, A i A i + h (u + *)AJAx 



i=2 



The real functions A, S, e are implicitly defined in terms of the auxiliary function rj(t) 



r/(t) = exp \ - J 



X(t')dt' -iwt-i / S(t')dt' , 



which satisfies the integrodifferential equation 



fj + iu)T] + 



/Vemv^'-^t) = o, 

JO u 



(9) 



(10) 



(11) 



subject to the initial condition rj(0) = 1. Moreover, considering the environment in thermal equilibrium, we have 



e(t) = 



|^)| 2 d /V MV(/3) 



2 dt 



E 



dre-^-^riir) 



(12) 



where nk(P) is the mean excitation number for the k th mode of the environment at inverse temperature (3 — 1/ksT. 
If the usual Born-Markov approximations hold, then S(t) = 0, X(t) = J^. ki := k, and e = kn, where ki characterize 
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the Markovian evolution when only the i th original oscillator is coupled to the bath, and n is the environment mean 
number of thermal excitations at frequency us. In this case the master equation simplifies to 



dp 



-£=-i"J2 [aJA^p] +k{n + 1) (2A 1 pA t 1 - AjA^ - pA^ 

i=l 

- kn ^2A|pAi — AiAjp — pAiAj 



(13) 



As one should expect by eq. (0), the master equation above describes the dissipative evolution of the collective mode 
Ai, and the independent unitary evolution of the remaining modes Aj, i > 2. It also should be noted that the 
damping constant k of mode Ai is larger (in general, much larger, for large N) than the individual constants fcj of 
the modes a^. This should be compared to the superradiance analogy discussed at section ITTTI 



V. TWO OSCILLATORS IN A DISSIPATIVE ENVIRONMENT 

From now on, we especialize to the case of two harmonic oscillators (N = 2). The collective modes can thus be 
written as 



AJ = cos(0)a} + sin(0)a|, AJ = - sin(0)a| + cos(0)aJj, 



(14) 



with tan0 = G%jG\ (this quotient can be taken as a positive number, if necessary, by redefining the mode a.2). If 
Markovian approximation is made, this relation takes the form tan0 = yk^Jki- The expression (1141) can be considered 
as giving |a|| by applying a rotation operator 71(0) on the set It is usual to call an operator which acts 

on operators a superoperator . Thus, 1Z (0) is a rotation superoperator . It is convenient to represent superoperators 
using the algebraic relations among the operators in which they act on, and a useful notation is to introduce a 

dot (•) in the position where the operator to be acted on must be placed. For example, 



A,. 



2 _ 
B 



B 



A, B 



A, 



A, B 



With this convention in mind one can verify that 



1l{6) = exp j(9 



aiaj, - a 2 a{, • 



and 



(15) 



Superoperators are very useful to study time evolution. As one can define an evolution operator U (t) by 
U (t) IV 7 (0)) = \i{s(t)), the evolution superoperator is defined by hi (t)p(0) = p(t). One completely solves the time 
evolution of a system by writing its evolution superoperator. The equation JSJ can be solved by the superoperator 



— iwt AIA 2 



3 (i-ti)A|»Ai e xA|Ai« e a:*»A|Ai f> zAi«A| 



U{t) = e -—r-a— <-ive^ i-±- e ~ — i—* e - 

where the coefficients v(t), x(t) and z(t) can be given in terms of the functions fj(t), eq. (|10f) . and Af(t), 



Af(t) = / dre(r) 
Jo 



(16) 



(17) 



as follows 



v(t) 



1 



l+Af(t) 



, x(t) = In • 



If the Markovian limit is applied, the preceding formulas reduce to 

^ e (-iw— fci -k 2 )t 



Z(t) = 1 



l + n(l-e- 2 ( fc i+ fe 2)*) : 



x = In ■ 



v /l + n(l-e- 2 ( fc i+ fc 2) t ) 



\v(t)\- 2 

1+Af(t)' 

(n + 1)(1 - e - 2 ( fc i+^)«) 
l + n(l-e- 2 ( fe i+ fe 2)*)) ' 



(18) 



The evolution superoperator U (t) can be expressed in terms of the original mode operators a^ by using the rotation 
superoperator 1Z (8) in the following way 



U(t) = e e [ a i a 2- 6 2a^»] e -iwt[a+a 2 ,»] ue (l-t>)a}«ai e za}ai» e ^^ 



(19) 
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Since the second collective mode is effectively decoupled from the environment, any density operator in the Hilbert 
space of this mode, times the asymptotic density operator of the coupled collective mode, provided it exists, will 
experience a unitary evolution. For simplicity we will further on restrict ourselves to the zero temperature case. Any 
density operator of the form (kets |m,n) refer to the original modes a^): 



E 



Pn 



E 



E 



(A|)"|0,0) (0,0| (a 2/ 

af sin 9 + k\ cos 9^j |0,0) (0, 0| (-&i sinfl + a 2 cos( 
Vnhn\(- sin e) n+m ~ ni ~ mi (cos 9) ni+mi 



(20) 



Pn,m 



n,m,ni ,mi 



(n — 7ii)!?ii!(m — mi)!mi! 
x \rii, n — n\) (mi, m — mi | , 

will be protected against dissipation and decoherence. In fact, applying the evolution superoperator to an initial 
density matrix of this form, we obtain 



E 



-iu)t(n-m) 



\fn\m\{— sin^ 



\n+m— m — mi 



(cos 9) 



n.m.ni ,mi 



w(ri — ni)!ni!(m — mi)!mi! 
x |ni,n — Tlx) {7711, m — mi| , 



(21) 



as one must expect. 

Now, we use the evolution superoperator on the initial operator density 

p(0) = (cos a 1 1,0) +e^sina|0, 1)) (cosa (1, 0| + e _i * sina (0, 1|) , 

which can be viewed as a one photon Fock state of the mode given by the creation operator 



(a, 4>) — cosa a} 



e lr ^ sina a^, 



(22) 



(23) 



where a and <p can be compared to Stokes parameters describing polarization. We will obtain how the dissipative 
properties of the mode (a, <f>) depend on these parameters. As this is a natural way for experimentally test DFS0, in 
the next section a more realistic situation is discussed. The above state asymptotically approaches the rank 2 density 
operator given by 



where the state 



p t _ >oc = P|^) M + (1-P)|0,0) (0,0| 
depends only on the individual decay rates ki, 

/fe|i,o)-Vfci"|o,i) 



and the weight P of this state is given by 



Ht2 cos a — \fk~i e 1 ^ sin a 



(24) 



(25) 



(26) 



Observe that varying a and 4> we can go from total preservation to total leakage. For example, if we set tan(a) = 
\fk~ijk~2, and = then the full state will leak to the ground state |0, 0), since in this case a = 9 and the initial 
photon was in the "superradiant" mode Ai. On the other hand, if we set tan(a) = —^Jk-z/ki, and </> = then the 
initial state will be exactly equal to \ip) (one photon in mode A2), and will persist at all times with probability 1 
(aside for an unimportant global phase). All other combinations will go to the density operator l|24() . which can be 
considered as an ensemble of pure state \ip) with probability P and the ground state |0, 0) with probability 1 — P. 
One can define the asymptotic fidelity, Poo (cc, 0), which is the overlap between the initial and asymptotic density 
matrices. In the above example it is given by 



Poo (a, 4>) = 



(^/k^cos(a) — \fk[e % ^ sin(a))(v / ^2 cos(a) — \fk[e ^sin(a)) 

ki + k 2 



(27) 
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VI. EFFECTS OF MORE REALISTIC MODELING 



We remark that the results above were obtained under a number of assumptions, which will be relaxed below. 
Notice that the use of the rotating wave approximation (RWA) is not essential in obtaining the decoupled mode: any 
interaction linear in the field operators would be as good (provided the other assumptions hold). Had we chosen an 
interaction linear in the identical oscillators but nonlinear on the environmental operators we would have obtained 
also a decoupled collective mode. In these cases, however, the complication would be only of technical nature leading 
to (much) more complex dynamics. 

Another important hypothesis to obtain DFS is that of identical frequencies of the original main oscillators. Of 
course, any interaction between them would destroy the symmetry upon which the existence of DFS rests. On 
the other hand, we have assumed that the oscillator-environment coupling satisfies gik — GiDk, which amounts to 
a separable coupling. It is not an easy task to find realizations of such interactions in nature given its nonlocal 
character. However, it might be a good approximation in special circumstances, as e.g. optical cavities. A particular 
consequence of the separability hypothesis can be seen writing the master equation, in the zero temperature limit, in 
terms of the original oscillators (with different frequencies for generality) 

Co = (—iu>i — ki) ajai • + (iu>i — k{) • ajai + 2/ciai • a{ + 
(— it02 — k 2 ) a^a 2 • + (iuj 2 — k 2 ) • a^a 2 + 2k 2 a 2 *a\ + 

k$ ^2ai • aj, — a^ai • — • ajjaij + k^* ^2a 2 • aj — a{a 2 • — • a|a 2 ^ . (28) 

The new quantity k^ appears since we consider the same environment interacting with both oscillators. These terms 
can be considered as an interaction between the oscillators mediated by the environment. If the separability condition 
is fulfilled, | /C3 1 2 = fcifc 2 , and if the oscillators are identical, eq. (|28|l is the same as eq. I|13[) . The other limit case is 
to consider both oscillators interacting independently with the environment. In this situation, the independence of 
the phases of the interaction coefficients gik will make their net effect on £3 null, and the eq. I|28l) will just describe 
two independent damped harmonic oscillators. Our interest is to study the above equation when the conditions for 
the existence of DFS are almost satisfied, i.e. {k^ 2 w k\k 2 and ui\ sa oj 2 . It must be noted that, by Cauchy-Schwarz 
inequality, I /C3 1 < k\k<x. 

We should poin t out that we are not using the usual approach of perturbation theory, of adding a small perturbation 
eH'. In Ref. |15|. the authors have shown that DFS are robust up to order e in the perturbation, and all orders in 
time, so our approach must be connected to second order (in e) perturbation. 

The explicit solution to this problem is given in the appendix. As is expected, there is no DFS without the 
separability and degeneracy assumptions, but if we are close to this conditions, we can obtain states much more robust 
to decoherence and dissipation than others. As in the previous section, consider the one photon states of eq. J22J>- 
Then we can define a weak decoherence mode (WD), which tends to the DFS when degeneracy and separability are 
approximated, and a strong decoherence mode (SD) which is analogous to the superradiant mode Aj. We want to 
explore the slight deviations from separability and degeneracy, so we define Sk and 5u> by 



Sk = yf~k\k 2 - N, 25u = u-2-ux, (29) 

and consider Su> <C w,: and Sk -C fcj. 

As in the previous section, varying the parameters of the initial state (1221) can be interpreted as varying the mode 
of the initial photon. In the regime discussed above we obtain: 

Ag D = J— (./MiJ ± ^e^Ta') , (30) 



wd \Jk\ J rk 2 \ V 2 

which must be compared to eq. (|23|l . The new parameter k is a kind of effective mean damping, and in the regime 
here discussed can be considered as k ss {k\ + k 2 ) /2. Each mode has its own damping constant, and this two are the 
extrema. Explictly we have 



for the weak decoherence mode, and 



2Sk^k 2 

kwD — — ; — : (31) 

k\ + k 2 



k S D = h+k 2 ^2k (32) 
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for the strong decoherence mode. One must note that while ksD is of the same order as the individual damping 
constants fej, the value of k\yD can be much lower. In the experiment with ionsQ it was exactly this lowering of the 
damping constant that was exhibited as an evidence of decoherence "free" subspaces. In the same experiment, one 
can see that the difference in damping constants is much larger in the situation with an engineered noise applied, 
since in this case the interaction with the environment is much closer to the separability condition. 

VII. CONCLUDING REMARKS 

We have studied decoherence and dissipation free modes for systems of harmonic oscillators. We discussed sufficient 
conditions for their existence. Although theoretically simple, these conditions are very difficult to be implemented in 
practical experiments. So we studied the slight deviations of this conditions, and instead of decoherence free subspaces, 
we obtained weak decoherence modes. This suggests that weak decoherence subspaces can be used to store quantum 
information for times much larger than individual carriers would be able to, even without being rigorous DFS. 

We compare DFS to the so called super and subradiance effects of a maser. In fact, when in the last section we 
compare the damping constants for weak and strong decoherence (eqs. (|31[) and (j22l) I this effect mimmimics inter- 
ference problem, where we are comparing the maximum and the minimum of a certain quantity in which interference 
effects arc recorded (in this case, the damping constant). 

It is important to stress that larger deviations of the rigorous conditions for DFS preclude the existence of even 
weakly decoherence subspaces, by making the decoherence time scales for such states smaller. However, one can 
conjecture that this kind of mechanism is so general that whenever an experiment obtain quantum mechanical results, 
it is testing some kind of DFS (e.g. the fullerenes experiment ) • 
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Appendix: Realistic model of two oscillators in details 

The evolution superoperator for ea. (|28[l can be expressed as pi) 



where 



0e" 



i(t)a^ai ■ m* (t)-ajai o(i)a+a 2 - q* (t)-af a 2 



R 



k 2 + fci i (lo 2 + ui ) 



c = k 2 ~ ki + i (lo 2 ~ , r = J~c 



c 2 +4:k§, A± = c±r 



q(t) = 2k 3 (l-e rt ) (A+e r 

,-R t 



for r^O 



= mi(t) _ 



2r 



-e-V (A+e r * - A_) , e" 1 *' = e" 2 * 



j 2 (t) = (l + \q(t)\ 2 ) 



n (t)= e- m ^+q(t) 2 e- m ^ + (\q(t)\ 2 ) 



,m 2 (t) 



- 1 



z(t) = -q(t)e- (m '^ +m ^ - q*(t) (f + |g(i)| 2 ) 



12 (t) 



(33) 

(34) 
(35) 
(36) 

(37) 

(38) 

(39) 



In this calculation neither the separability nor the degeneracy (even approximated) conditions have been used so far. 
For the sake of comparison we use the same initial condition of section (eq. (|22J) V In the general case its time 
evolution is also given by 



pit) = p(t) \m) m)\ + (i - p(t)) io, o) (o, oi 



(40) 
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but now the state \ip(t)) is, 

, _ (co S (g)M_(t) + S in(g)e^Q(t))|l,0) + (sin(g)e^M + (t)+cos(g)Q(t))|0,l) 

■ (41) 

and its coefficient is given by 

P(t) = |cos(0)M_(t) + sin(6>)e^Q(t)| 2 + |sin(f)e^M + (t) + cos(6>)Q(t)| 2 , (42) 
where the functions M±(t) and Q(t) are given by 

M±(t) = ^ (e- r */ 2 (l T -) + e rf / 2 (l ± -)) , Q(t) - ^e~ flt (e^ 2 - e rt ' 2 ) . (43) 

Now, we assume slight deviations from degeneracy and separability, that is, u)\ = uj — Suj, lu 2 = lo+Slo, k% = y/kikz — Sk, 
with Su <C ui, 5k <C \Jk]k 2 . Then, the state \ip(t)) can be approximated as 

m) = e -- ( 6MWi e -w + 6|i,o) + 6|o,i) e -2^A (44) 

where Q , £j do not have any temporal dependence and are given by 



{ki =F iSiv) cos(a) ± \fk\k 2 sin(a)e 1 ^ 

^ = — rrr > ( 45 ) 

2 ki + k 2 

(k 2 ± i^w)e l< ^sin(a) ± \fk\k 2 cos(a) 
& = — 7—— r ■ (46) 

2 Kl + «2 

In the general case it is not possible to find initial conditions which are completely decoherence free. Nor it is possible 
to find two orthogonal subspaces with very different characters in what decoherence is concerned. However, we 
can choose the initial condition as to have a minimal component either in a strong decoherence (SD) or in a weak 
decoherence (WD) subspaces, by choosing, e.g. 



tan{a) sd = ± A 



k. 



^ ki wd k 



where k is some average dissipation constant. The corresponding states, apart from a phase, can be written as 



VSD ) = 1 ■ ( .fh |1,0) ± Ik^^ |0,1)) 

WE> / Vki+k 2 \V 2 V l / 



(48) 

If (5w <C ki,k 2 the phase can be ignored. Moreover, if k\ = k 2 then we have k = fci = &2 and the phase can 
be unaunambiguously determined. The weak decoherence wavefunction defines a mode which is robust against 
decoherence. The damping constant of this mode can be read from eq. (|44|) as the value given in eq. I|31|) . Analogously 
for the strong decoherence mode, with damping constant given by eq. 132fl . 
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